Chiral and deconfinement phase transitions of two-flavour QCD 
at finite temperature and chemical potential 
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We present results for the chiral and deconfinement transition of two flavor QCD at finite tem- 
perature and chemical potential. To this end we study the quark condensate and its dual, the 
dressed Polyakov loop, with functional methods using a set of Dyson-Schwinger equations. The 
quark-propagator is determined self-consistently within a truncation scheme including temperature 
and in-medium effects of the gluon propagator. For the chiral transition we find a crossover turning 
into a first order transition at a critical endpoint at large quark chemical potential, Hep/Tep ^ 3. 
For the deconfinement transition we find a pseudo-critical temperature above the chiral transition 
in the crossover region but coinciding transition temperatures close to the critical endpoint. 



Introduction 

Heavy ion collision experiments at RHIC, the LHC 
and the future FAIR project are designed to study 
the properties of the large temperature phase of 
QCD, the quark-gluon plasma. Of particular in- 
terest are the details of the chiral and deconfine- 
ment transition. Here, a standard scenario favors 
a chiral crossover at small chemical potential turn- 
ing into a first order chiral transition at a critical 
point. The search for this critical point is one of 
the main motivations in the physics program of the 
future CBM/FAIR experiment in Darmstadt and 
the NICA project in Dubna. However, on the the- 
oretical side there is ambiguity not only about the 
location of the critical point but also if the stan- 
dard scenario is even correct. Among other possi- 
bilities [ij an exotic, quarkyonic matter phase [2| or 
inhomogeneous chiral condensates 0, Q could be 
realized. Unfortunately, lattice Monte-Carlo simu- 
lations cannot be used to resolve this issue due to 
the notorious sign problem. 

Effective models to QCD like the Polyakov- 
loop extended Nambu-Jona-Lasinio model (PNJL) 
and the Polyakov-loop extended quark-meson 
model (PQM) 0, 0] constitute an alternative ap- 
proach. Much progress has been made taking into 
account fluctuations beyond mean- field 0, which 
shift for example the critical endpoint to larger 
chemical potential. Interestingly this is in agree- 
ment with recent lattice results indicating that the 
critical endpoint if existent is at much larger chem- 
ical potential [§| than previously thought Q. 

Within QCD fortunately also non-perturbative 
tools not limited by the sign problem are at our 
disposal. In particular, it has been shown that 
Dyson-Schwinger equations (DSEs) and the func- 



tional renormalization group are capable 

to describe both, the chiral and the deconfine- 
ment transition at finite temperatures and imag- 



inary chemical potential |12l. Il4|. 



In this letter we extend these applications to fi- 
nite, real chemical potential. We use the DSEs for 
the propagators of Landau gauge QCD with two 
fermion flavors in a suitable truncation scheme. 
Working with physical quark masses we determine 
the quark condensate and its dual, the dressed 
Polyakov loop in the real (T, ^)-plane. We ex- 
tract transition temperatures for the chiral and 
deconfinement transition and determine their na- 
ture. We locate the critical endpoint of the chi- 
ral transition. Our study provides the tools for 
further, more detailed investigations of gauge in- 
variant properties of the QCD phase diagram by 
functional methods. 

Order parameters 

We calculate the order parameters for chiral sym- 
metry breaking and confinement from the fully 
dressed quark propagator, given by 

S-^{p) ^ i^pA{p) + i-n{LOr, + i^l)C{p) + B{p), (1) 

where is the quark chemical potential and cj„ = 
7rr(2n + 1) are the Matsubara modes in the imagi- 
nary time formalism with temperature T . The ar- 
gument (p) serves as an abbreviation for (]y^,ujn)- 
The scalar and vector dressing functions B and 
A, C are calculated from the quark DSE. For the 
bare, renormalized propagator Sq^{p) we have 
A ^ C ^ Z2 and B = m, where m is the renormal- 
ized bare quark mass and Z2 is the quark renor- 
malization factor. 

There are several equivalent choices of order pa- 
rameters for chiral symmetry breaking; here we use 
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the quark condensate 



(#) ^Z.Tj^l (1^3 Trz, [Sip)] . (2) 



We work with two degenerate quark flavors with 
physical quark masses and therefore expect a 
smooth crossover at vanishing chemical potential 
0. Below, we determine the transition tempera- 
tures via the maximum of the susceptibility ^^^^ . 

Extracting an order parameter for confinement 
from the quark propagator is a much more chal- 
lenging task. In the last years a method has been 
developed to calculate dual condensates I]„ which 
are sensitive to centre symmetry breaking 17 - 1^ . 
They are defined by the transform 



where = J^^ ^ and (ip'>p)cp is the quark conden- 
sate evaluated at generalized, ?7(l)-valued bound- 
ary conditions tp{x,l/T) = e*'^^/'(x, 0) with ip G 
[0, 2tt[. On a lattice, the quantity S„ can be in- 
terpreted as a sum of all closed loops winding n 
times around the Euclidean time direction. Con- 
sequently, S+i contains the Polyakov loop and S-i 
its conjugate, together with all other oriented loops 
that wind once around the space-time torus. They 
therefore have been called 'dressed Polyakov loop' 
and its conjugate [l^l- For to act as order 
parameters for deconfinement it is mandatory to 
implement the generalized, ?7(l)-valued boundary 
conditions only on the level of observables, but not 
in the partition function itself. All closed quark 
loops therefore maintain the physical value </5 = tt, 
whereas ip G [0,27r[ for the test quark in Eq. (jS]). 
This procedure breaks the Roberge- Weiss symme- 
try, a necessary condition for the dual condensate 
to act as an order parameter for centre symmetry 



(3) 



breaking [IJ, [2 



For our choice of physical quark mass we ex- 
pect a deconfinement crossover at small /i. To 
determine the pseudo-critical temperature we use 
the maximum of the derivative with respect to the 
quark mass. 

Dyson-Schwinger equations 
The quark propagator is calculated from its Dyson- 
Schwingcr equation (DSE), shown in the first line 
of Fig. [1] It depends on the dressed quark-gluon 
vertex and the dressed gluon propagator, which 
need to be specified. For the Landau gauge gluon 
propagator at finite temperature, the most reli- 
able source up to now are quenched lattice calcula- 
tions Il3l l21i . Since we are interested in studying 




FIG. 1: The DSE for the quark and gluon propagators. 
Filled circles denote dressed propagators and vertices. 
The shaded circle denotes the quenched ghion. 



transitions at non-vanishing chemical potential it 
is mandatory to include quark in-medium effects in 
the gluon propagator. To this end wc use the lat- 
tice data from Ref. [31 for the quenched gluon and 
add the quark polarization tensor from the gluon 
DSE to effectively unquench the system, see sec- 
ond line of Fig. [T] This procedure is not entirely 
self consistent. In particular unquenching effects in 
the Yang-Mills part of the gluon-DSE, i.e. in the 
ghost-loop and gluon-loop diagrams hidden in the 
quenched gluon in Fig. [T] are neglected. In the 
vacuum wc compared these effects with the fully 
selfconsistent treatment of Ref. [g^]. We found 
that the approximation used here overestimates 
the unquenching effects, with momentum depen- 
dent deviations in the gluon propagator below the 
five percent level. Assuming that this remains so 
at finite temperature and chemical potential one 
would expect to slightly underestimate the transi- 
tion temperatures perhaps on the 5-10 MeV level. 
We believe that such a treatment is sufficiently ac- 
curate for the qualitative study envisaged in this 
work. 

Within the polarization tensor we additionally 
neglect the quark self-energies in analogy with cor- 
responding works on finite density- Q CD in the 
DSE formalism [i^. This approximation is justi- 
fied by findings comparing lattice and DSE results 
for the quark spectral function 15[ showing that 
the self-energies are suppressed close to the transi- 
tion temperature. The tensor is then given by 



^ Tr [5o(g)7p5o(A;)rO] 



(4) 



(27r)3 



the vertex 



with g = P + fc, ^fc = E„/ 

renormalization Z\p and the quark-gluon vertex 
for bare propagators r°(p) = g^yV(j>), discussed 
below. 

Projected to modes longitudinal (L) and 
transversal (T) to the medium, the gluon DSE 
reads 



the transition temperatures and the nature of the 
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with the quenched and unquenched gluon dress- 
ing functions Z!Pj^ and Zt,l, respectively, and the 
quark loop 117- ^. Equation ([S]) may be evaluated 
in leading order hard thermal loop where we obtain 
thermal momentum dependent masses 



NfTiaT^hip) 



with arxip) = ZiFg^ZT,L[p)T°{p)/{ZlA^). 



For the quark-gluon vertex at finite temper- 
ature no information is available yet, neither from 
the lattice nor from functional methods. For the 
purpose of this study we therefore rely on a phe- 
nomenological ansatz that has been developed in 
earlier works, Ref. 1^, llSf. To make this work 



self-contained we repeat the vertex here, 



C{k) + C{p) 




r 



/3oa(/i)ln[gVA2 + l] 



2<5^ 



47r 



(7) 



where q ~ {q,ujq) denotes the gluon momentum 
and p = (p,ujp), k = {k,ujk) the quark and anti- 
quark momenta, respectively. Furthermore 25 = 
— l8Nc/ {iANc — SNf) is the anomalous dimension 
of the vertex. This ansatz contains the quark dress- 
ing functions A and C in a manner dictated by the 
Slavnov- Taylor identity of the vertex. It therefore 
implicitly depends on temperature and chemical 
potential in a meaningful way and also contains 
unquenching effects via the corresponding changes 
in A and C. In the quark-DSE we keep this de- 
pendence, whereas in the gluon polarization we set 
^4 = C = 1 to be consistent with the HTL approx- 
imation described above. 

The parameters di^2 and A of the vertex have 
been fixed in Ref. [13| for the case of quenched 
QCD. The temperature independent scale A is as- 
sociated with the scale of the lattice data for the 
gluon propagator and given by A ^ 1.4 GeV. 
The parameter d2 is associated with the transi- 
tion from the logarithmic ultraviolet running of 
the vertex towards the infrared constant and is 
given by ^2 = 0.5 GeV^ Moreover, di = 4.6 GeV^ 
parametrizes the infrared strength of the vertex. 
Both, di and ^2 have been fixed in Ref. [l^l such 
that the chiral and deconfinement transition tem- 
peratures extracted from the quenched quark prop- 
agator matched the deconfinement transition tem- 
perature taken from the Yang-Mills lattice input. 
Note that this procedure is not very sensitive, i.e. 
it does not require any fine-tuning: once di and c?2 
are in the right ballpark, the resulting transition 
temperatures have been found to be insensitive to 
variations of di and d2 up to the ten percent level. 
Thus once the vertex is chosen such that it roughly 
matches the lattice gluon propagator, its details 



arc not so relevant. This is certainly reassuring. 
Using a bare quark mass of m(/i = 30 GeV) = 3.7 
MeV we also determined the resulting (quenched) 
pion mass at zero temperature, which is roughly 
ten percent above its physical value and indepen- 
dent of changes in di and c?2- 

In our unquenched calculation we use the same 
vertex parameters as in the quenched case. This 
choice is only justified by simplicity. In general, 
one also expects unquenching effects in the vertex, 
which may affect the transition temperatures. One 
way to determine the corresponding changes in di 
and d2 would be to calculate unquenched observ- 
able like the pion mass and decay constant in the 
T = limit and match to experiment. Unfortu- 
nately this cannot be done within the HTL approx- 
imation of the quark polarization used in this work, 
since it only contains temperature dependent fluc- 
tuations which vanish in the T = limit. A more 
refined treatment of the quark backreaction includ- 
ing also quantum fluctuations is under way and will 
be presented in a future work. We do, however, not 
expect that unquenching effects in the vertex will 
have a large impact on the transition temperatures. 
At zero temperature and with fixed gluon propa- 
gator these effects have been found to reduce light 
hadron masses on the ten percent level [24|. Thus 
with fixed gluon propagator one may expect that 
unquenching effects in the vertex will reduce the 
transition temperatures by a similar amount. This 
is, however, counterbalanced by a corresponding 
decrease of quark polarization effects in the gluon 
propagator, which in turn increases the transition 
temperatures again. Thus the combined effect may 
indeed be very small. 
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FIG. 2: The quark condensate as well as the dressed 
Polyakov-loop as a function of temperature in two- 
flavour QCD. 



FIG. 3: The quark condensate {'4i'4>)v ^ ^ function 
of temperature and boundary angle ip at fixed quark 
chemical potential /x = 100 MeV. 



Results 

As a first check of the employed truncation scheme 
we determine the transition temperatures aX ji — Q. 
In the quenched system, i.e. without the quark- 
loop contribution of Eq. our system of equa- 
tions reduces to the one outlined in Ref. and 
we obtain Tc'=° = T^J^^^j = 277 MeV. This 
is in accordance with lattice caluclations. When 
we switch from Nf = to Nf = 2 the main ef- 
fect of the quark loop onto the gluon propagator 
is to decrease its strength and to lower the char- 
acteristic scale Aqcd- This also results in de- 
creased chiral and deconfinement transition tem- 
peratures. We obtain T^'"^ = 180 ± 5MeV and 
'^deconf ~ ^^"^ ^ SMcV which is again close to the 
corresponding values on the lattice [IHj. Here the 
error bar reflects possible numerical uncertainties 
in our calculation. The additional sources for sys- 
tematic errors of our calculation have been dis- 
cussed in the last section. As a rough guidance 
for the reader we guesstimate an additional sys- 
tematic error of ±20MeV. 

Our results for the quark condensate and the 
dressed Polyakov-loop are shown in Fig. [2l In the 
quenched calculation of Ref. [l^l a sharp rise and 
fall of the quark condensate close to Tc has been 
observed, which has been attributed to a corre- 
sponding behavior of the electric part of the gluon 
propagator. Here, the backreaction of the quark- 
loop has suppressed this effect and the conden- 
sate decreases monotonically as expected. This 
change of behavior together with the correct order 
of the transition temperatures gives us confidence, 
that our implementation of the backreaction of the 
quarks onto the Yang-Mills sector serves its pur- 
pose. 

We then switch on the chemical potential and 
determine the chiral condensate as well as the 
dressed and the conjugate dressed Polyakov loop 



in the (T, /i)-plane. Note that for /i ^ 0, {•ipip)ip de- 
velopes an imaginary part, which will be discussed 
below. A typical example for the behavior of the 
real part of the (/3-dependent quark condensate at 
a moderate value of fi is shown in Fig.[3l At = tt 
we observe the chiral cross-over with temperature 
indicated by the ordinary quark condensate. At 
if = 0, 27r the condensate rises with temperature 
similar to = as discussed in Ref. [12|. The 
deconfinement transition is observed along the ip- 
direction, when the constant behavior of (V'V')v 
below Tc changes into the typical variations ob- 



served above Tc jl2| . The new element at finite fi 



is the behavior for even larger values of T: There 
{ipip)ip develops an additional discontinuous struc- 
ture along the (/^-direction. We have checked that 
this structure is genuine wrt. variations of our ver- 
tex ansatz. For increased chemical potential the 
temperature where this discontinuity opens moves 
closer to the deconfinement transition tempera- 
ture. This suggests that the deconfinement tran- 
sition may turn from crossover to second or first 
order at chemical potentials as large as the one of 
the chiral critical endpoint (see below). Within the 
accuracy of our numerical calculations we could, 
however, not yet resolve this issue in detail. This 
will be addressed in future work. The imaginary 
part of not shown in the figure, is antisym- 

metric wrt. to </3 — > 277 — (/J. Therefore Ej-i = 

{Re{-ip'tp)^ cos{ip) ± Im^tpip)^ sin((y9)} stays real 
and the difference between quark and anti-quark 
loops is generated by Im{ipil')v 

In Fig. UJ we display the main result of this 
work: the phase diagram of two-fiavour QCD. 
We find a chiral crossover at small chemical po- 
tential. The deconfinement crossover happens at 
somewhat larger temperatures and one observes a 
characteristic (albeit small) splitting of the tran- 
sition temperatures from the Polyakov-loop and 
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FIG. 4: The phase diagram for chiral symmetry break- 
ing (x) and deconfinement of quarks (Ei) and anti- 
quarks (E_i). 



its conjugate At larger chemical potential 

all transitions come together again. The chiral 
crossover line goes over into a critical point at ap- 
proximately {Tep, fJ.Ep) ~ (95, 280) MeV, followed 
by the coexistence region of a first order transi- 
tion. Thus we find the comparatively large value 
I^ep/Tep ~ 3. We have checked, that these val- 
ues are not overly sensitive to the details of our 
truncation: when changing the parameters in the 
vertex ansatz within a reasonable range we observe 
variations of {Tep, (J-ep) of the order of ten per- 
cent. Thus a firm conclusion of the present ap- 
proach seems to be that hep /Tep ^ f: If there 
is a critical endpoint, it happens at large chemical 
potential. This statement agrees with the result 
of corresponding calculations in the PQM model, 
once quantum corrections have been taken into ac- 
count Q- Expectations from recent lattice calcu- 
lations at Nf = 2 + 1 close to the continuum limit 
also seem to point in this direction Q . 

Certainly, at such large values of the chemical 
potential, our truncation scheme may no longer 
be reliable: baryon effects that are not implicitly 
included in our truncation of the quark-gluon in- 
teraction may play an important role here. Also 
the formation of inhomogcneous chiral condensates 
may be favored upon the homogeneous one studied 
here. We believe that our work provides a suitable 
basis for further investigations in these directions. 
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